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Se presenta el diseño de un sistema de diagnóstico basado 
en relaciones redundantes y observadores de estado 
—que funciona manteniendo al sistema de transporte 
en oOperación— para la detección, identificación y 
reconstrucción de fallas en ductos, tales como: anomalías 
en los instrumentos de medición, bombas y extracciones 
desconocidas. Los algoritmos se desarrollaron a partir 
de ecuaciones diferenciales parciales no lineales que 
caracterizan el comportamiento del fluido —considerando 
el principio de conservación de masa y momento— y que 
fueron discretizadas en el espacio utilizando el Método de 
Diferencias Finitas. Con la finalidad de discernir entre los 
diferentes tipos de fallas y de reconstruir el comportamiento 
de éstas, el sistema de diagnóstico opera por etapas. La 
primera etapa, llamada de Detección y Aislamiento de Fallas, 
tiene como objetivo aislar el síntoma de la falla con ayuda de 
un conjunto de relaciones redundantes deducidas a partir 
del modelo nominal del ducto. Esta etapa de aislamiento 
simplifica la segunda etapa: Reconstrucción de la Falla, en la 
que algoritmos de observación estiman la evolución temporal 
de la falla aislada. El sistema completo de diagnóstico se 
validó mediante una serie de experimentos realizados en 
un modelo hidráulico piloto de aproximadamente 200 m 
de largo, instrumentado y construido específicamente para 
llevar a cabo pruebas de monitoreo de ductos en el Instituto 
de Ingeniería de la Universidad Nacional Autónoma de 
México. 
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Introducción 


La supervisión y el monitoreo automático de 
ductos es un reto para la ingeniería debido a 
las consecuencias catastróficas que pueden 
desencadenar los ductos en mal estado, como son 
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This paper presents the design of an online diagnostic system 
(functioning while the transportation system remains in operation) 
to detect, identify and rebuild pipeline faults based on redundant 
relations and state observers. The faults include those occurring in 
measurement instruments, pumps, and unknown extractions from 
the pipeline. The algorithms that compose the diagnostic system 
were developed from nonlinear partial differential equations that 
characterize the behavior of the fluid according to the principle 
of conservation of mass and momentum. These equations were 
approximated in space using the Finite Differences Method. In order 
to distinguish between the different types of faults and reconstruct 
their behaviors, the diagnostic system operates in stages. The first 
one —which is called Fault Detection and Isolation— aims to isolate 
each fault symptom with the help of a set of redundant relations 
that are deduced from the nominal model of the pipeline, i.e., in 
normal conditions. This first stage simplifies the second, called Fault 
Reconstruction, which is composed of observation algorithms that 
estimate the temporal evolution of the isolated faults. The overall 
diagnostic system is validated through a series of experiments 
carried out in a pilot hydraulic pipeline at the Engineering Institute, 
UNAM. This pipeline is approximately 200 m long and was 
implemented and built specifically to carry out tests to monitor 
pipelines. 


Keywords: Diagnostic systems, fault detection, hydraulic pipelines, 
pipeline monitoring. 


la pérdida de vidas humanas y daños ambientales 
con alto impacto ecológico. En general, el objetivo 
de un sistema automatizado de supervisión 
en tiempo real consiste en localizar, tan rápido 
como sea posible, la presencia de fallas con un 
mínimo de instrumentación y costo. 
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A partir de sistemas SCADA (acrónimo de 
Supervisory Control and Data Acquisition) 
que operan en tiempo real, donde se registran 
presiones y flujos en puntos limitados de las re- 
des de ductos, se han implementado localiza- 
dores automáticos de averías utilizando 
sensores físicos y virtuales, diseñados a partir 
de algoritmos de estimación e identificación 
(Stoianov, Lama, Sam, Timur, € Csail, 2007; 
Korbicz, Koscielny, Kowalczuk, € Cholewa, 
2004; Makar £: Chagnon, 1999). En los ductos, 
las anomalías de mayor interés a localizar 
en tiempo real son fallas en bombas y servo- 
válvulas, fugas en el cuerpo del ducto y 
obstrucciones. Si los ductos están bajo tierra o 
agua, la inspección se dificulta, pero además, 
cuando no se puede suspender la operación del 
ducto, el diagnóstico debe hacerse considerando 
los fenómenos dinámicos asociados con el 
fluido en cuestión. La detección de fallas de 
magnitud considerable resulta más simple 
que la de pequeñas fallas, pues los sensores no 
siempre tienen la sensibilidad requerida para 
indicar desviaciones pequeñas, además de que 
el ruido presente en las mediciones afecta la 
precisión de los estimados. 

En cuanto al diseño de los algoritmos 
para la detección y localización de una fuga 
en ductos se han usado diversos modelos 
de comportamiento estático (Billmann «: 
Isermann, 1987; Carrera €: Verde, 2001). De 
manera reciente, estos modelos se ha extendido 
con éxito al caso de ductos con perfiles 
topográficos no uniformes (González, Verde 
y Torres, 2013). No obstante su simplicidad, 
la principal desventaja de los algoritmos 
basados en estos modelos es que no pueden 
utilizarse si el fluido no opera en condiciones 
aproximadamente estáticas. 

Como solución alternativa, Brunone y 
Ferrante (2001) formularon una metodología 
a partir del comportamiento de la onda 
de presión que se refleja hacia su punto 
de origen al encontrar una fuga. En estas 
condiciones, el tiempo de arribo de la onda 
permite determinar la posición de fuga. 
Wang, Lambert, Simpson, € Vtkovsky (2005) 


también estudiaron las características de la 
onda de presión y reportaron que se atenúa en 
función de la posición de la fuga. Por otro lado, 
Mpesha, Chaudry, $: Gassman (2001), y Covas, 
Ramos, € De Almeida (2005) desarrollaron 
métodos basados en la respuesta frecuencial de 
la tubería, cuya respuesta dinámica se puede 
obtener provocando transitorios periódicos 
mediante la apertura y el cierre de una válvula. 
Si hay fugas en el ducto, la respuesta frecuencial 
resonantes 
adicionales; a partir de la frecuencia de estos 


de la presión muestra picos 


picos puede estimarse la posición de la fuga. 
Ferrante y Brunone (2003) también encontraron 
un enfoque en el dominio de la frecuencia 
y resolvieron las ecuaciones del Golpe de 
Ariete utilizando el método de la respuesta al 
impulso. Mediante un análisis armónico de la 
solución, concluyeron que a partir de la presión 
transitoria se pueden identificar la posición y el 
tamaño de la fuga. 

Cuando se dispone de modelos analíticos 
precisos es posible obtener síntomas de fallas 
de sistemas complejos por software de forma 
muy rentable, usando las herramientas desa- 
rrolladas por la comunidad de control auto- 
mático (Isermamn, 2011). La clave para generar 
síntomas de las fallas es identificar incohe- 
rencias en el sistema a través de mediciones 
y evaluando un residuo, cuyo valor es cero 
en condiciones normales y se desvía de éste 
ante eventos anormales. Dicha incoherencia se 
puede obtener a partir de modelos matemáti- 
cos con la ayuda de filtros o relaciones de pa- 
ridad (Chow éz Willsky, 1984), estimación de 
parámetros (Isermann € Muenchhof, 2011), y 
técnicas de observación (Chen éx Patton, 1999). 
Estas herramientas explotan de forma directa 
O indirecta el comportamiento del sistema, al 
sustituir la multiplicidad de sensores físicos 
por algoritmos que se ejecutan en tiempo real. 
Para profundizar sobre el tema de EDI, acróni- 
mo de Fault Detection and Isolation, consultar 
Verde, Gentil y Morales (2012). En particular, la 
detección de fallas en la instrumentación de un 
proceso se ha atacado ampliamente con la idea 
propuesta por Frank (1990) de usar un banco 
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de observadores con tareas específicas. Recien- 
temente se presentaron algunas técnicas gene- 
rales para estimar las fallas en Zhang (2011) y 
Boulkroune, Galvez-Carrillo y Kinnaert, (2011). 

Los hechos arriba descritos motivaron 
el desarrollo del presente trabajo, en el que 
la principal contribución es un sistema de 
diagnóstico que considera de manera global 
fallas en instrumentos de medición, averías 
en bombas y fugas en el ducto. El sistema usa 
un algoritmo con capacidad para distinguir 
entre las diversas fallas en tiempo real con 
cinco residuos y un banco de observadores 
dinámicos de altas ganancias para reconstruir 
la evolución de la falla detectada. 


Preliminares de diagnóstico de fallas 


En la comunidad de FDI, una falla se define 
como una desviación no permitida de al menos 
una propiedad característica de un sistema 
(Isermann, 2006). Estas desviaciones pueden 
ocurrir en sensores, en actuadores o en los 
componentes del proceso. Al procedimiento de 
supervisión en tiempo real de un proceso con 
capacidad para detectar, aislar fallas y evaluar 
su significado y gravedad se le conoce como 
sistema de diagnóstico de fallas. Estos sistemas 
involucran tres tareas (Patton, Frank, e Clarke, 
2000). 


e Detección de la falla: tomar una decisión 
binaria sobre el estado del sistema de 
acuerdo con las condiciones nominales de 
operación, por ejemplo, establecer si en el 
proceso está presente un evento anormal. 

e Aislamiento de las fallas: determinar el tipo 
de falla presente, para lo cual se debe contar 
con un síntoma, como estipular qué sensor, 
actuador o componente se encuentra en 
condición anormal. 

e Identificación o reconstrucción de la 

falla: estimar el tamaño y tipo de la falla, 

incluyendo su evolución temporal. 


La idea básica para la detección es el uso 
de modelos de conocimiento con redundancia, 


para verificar que el comportamiento actual del 
sistema es consistente con el modelo nominal 
establecido. Una de las ventajas de disponer 
de redundancia de información es que ésta 
se puede evaluar simplemente mediante 
un proceso, bajo condiciones de operación 
sin necesidad de 
instrumentación física adicional. Sin embargo, 


bien caracterizadas y 


se debe tener un modelo de comportamiento 
adecuado del sistema. La figura 1 describe 
el principio de la detección, comparando 
modelos e identificando la incoherencia como 
consecuencia de la falla. Para generar redun- 
dancia analítica, el modelo dinámico de 
parámetros concentrados más común es la 


representación en espacio de estados: 


x(t)=f(x(t),u(t)0(0) 


y(0)=h(x(0) (0 
donde x(t)ER” son los estados del sistema; 
u(t)ER”, las entradas exógenas o señales 
de control; y(t)ER”, las salidas y variables 
medibles del proceso, mientras que 0(t)ER”o 
representa los parámetros físicos del sistema 
(como longitud, elasticidad, fricción, etcétera). 

Suponiendo conocido el modelo (1) y el 
vector de entradas u(t), se puede calcular la 
salida /(t); y si además, se dispone de la salida 
a través de mediciones y(t), se tiene dos veces 
la misma información y hay redundancia. En 
condiciones anormales del modelo, la salida 
medida y(t) es incoherente con la calculada y(t) 
y la redundancia se pierde. A las ecuaciones 
que se pueden obtener con el modelo y las 
mediciones con redundancia se les conoce 
como relaciones redundantes analíticas y éstas 
permiten generar residuos r(t) con la siguiente 
propiedad (Gertler, 1991): 


r(1)= O condición normal 0) 


0 condición anormal 


Así, r(t) diferente de cero denota error en el 
modelo matemático y con la ayuda de dichos 
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Señales conocidas 


Modelo de 
comportamiento 
normal 


Estado 
calculado 


Evaluación 
y decisión 


Síntoma 


Figura 1. Principio de la obtención del residuo r(t) y síntoma de fallas donde K(t) son las fallas de interés 


y K(t) denota perturbaciones y ruido. 


residuos se determinan síntomas para cada 
falla vía funciones de decisión. La figura 1 
describe este principio de detección usado para 
generar residuos y síntomas. 

En la mayoría de las aplicaciones no existe 
un residuo r(t) sensible solamente a una falla, y 
robusto al resto de ellas y a las incertidumbres 
en el modelo. En esta condición, es necesario 
obtener patrones de residuos para aislar las 
fallas a partir de análisis específicos de las 
ecuaciones del modelo y su comportamiento 
ante las fallas. 

La etapa de reconstrucción de fallas, una 
vez que se ha determinado cuál está presente, 
se basa fundamentalmente en algoritmos de 
estimación de parámetros y de observación 
de los estados. Las técnicas de estimación de 
parámetros son apropiadas si las fallas están 
asociadas con cambios en los parámetros del 
modelo. Es decir, se hace uso del hecho de que 
las fallas del sistema dinámico se reflejan en 


parámetros físicos. La idea de este enfoque es 
reconstruir las fallas por medio de la estimación 
de parámetros del modelo matemático (Ljung, 
1999). 

La técnica de observadores es apropiada 
si las fallas producen cambios en actuadores, 
sensores O variables de estado no medibles. 
Recientemente, se ha propuesto reconstruir las 
fallas usando un observador del modelo (1) 
en condición anormal, en donde el vector x(t) 
se extiende con los parámetros asociados con 
las fallas. Es decir, se considera que el vector 
extendido contiene los parámetros asociados 
con las fallas O(t)€ER"o formando: 


(3) 


donde x (t)€ER”*"w denota el espacio extendido. 
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Una propiedad de este tipo de estimación 
es que permite estimar conjuntamente 
parámetros y estados. Esta ventaja llevó a 
seleccionar a los algoritmos de observación 
para la reconstrucción de fallas en el caso de 
ducto. 


Algoritmo de observación 


Definición 1. Un observador de estados es un 
algoritmo dinámico, cuyo propósito es estimar las 
variables de un proceso x(t) utilizando el modelo 
(1), las mediciones disponibles (u(t) y y(t)), y un 
término de corrección del error entre las mediciones 
y los estimados, que garantice la convergencia de los 
estados. Dado que los observadores estiman estados 
por software, se les denota también como sensores 
virtuales (Besancon, 2007). 

El diagrama de bloques de un observador 
conectado a un proceso tiene la estructura 
mostrada en la figura 2, donde al error de 
estimación, e(t) = y(t) — y(t) corresponde a la 
diferencia entre la salida medida del proceso 
y(t) y la salida estimada (t). El objetivo del 
observador es hacer que el error e(t) converja 
a cero cuando el tiempo tiende al infinito. 
En la práctica, se desea que el error sea lo 
suficientemente pequeño en un corto tiempo. 


Una contribución importante en la teoría 
de observadores no lineales se reporta en 
Gauthier, Hammouri y Othman (1992). En 
dicho trabajo se presenta el diseño de un 
observador, nombrado de alta ganancia, para 
sistemas no lineales compuestos de un término 
dinámico lineal y de un término no lineal con 
forma triangular, acotado y afectado por las 
entradas (Khalil, 2002). La sintonización del 
observador, en este caso, se realiza ajustando un 
solo parámetro A, mientras que la convergencia 
está asegurada mediante la solución única 
de una ecuación de Lyapunov (Bornard éz 
Hammouri, 1991). 

El Apéndice A describe el algoritmo genérico 
del observador de alta ganancia propuesto por 
Gauthier et al., y que se emplea con modelos 
del fluido y de las fallas para reconstruir los 
parámetros de los escenarios de fallas del ducto 
y de la instrumentación. 


Modelo dinámico del fluido 


Asumiendo que la distribución de velocidad es 
uniforme, que el flujo es unidimensional, y que 
la densidad del líquido y área son constantes, la 
dinámica de un fluido en un ducto horizontal 
cerrado, sin tomas laterales, está gobernada por 
las ecuaciones de movimiento y continuidad: 


Estados estimados xt) 


Proceso 
Entradas 


Salida estimada 00) 


Salida del proceso 


Figura 2. Esquema de estimación con un observador de estados. 


Tecnología y 


Ciencias del Agua, vol. V, núm. 4, julio-agosto de 2014 


Torres et al., Algoritmos de diagnóstico para fallas en ductos 


nz) y? 9Q(z,t) SÍ (4) 
ot gÁA,  0z 


F 


1 0Q(21), oH(21) fQ(2,+)0(2,+) 
A, 0 Saz 244? 


r 


=0 (5) 


con (z,t) € (0,L) x (0, oo) como coordenadas del 
tiempo (s) y el espacio (m). Los parámetros que 
intervienen en el conjunto (4) y (5) se enlistan 
en el cuadro 1. 

Dado que no existe una solución analítica 
del sistema (4) y (5), es necesario utilizar 
métodos numéricos para obtener una solución 
aproximada. Algunas de estas soluciones 
pueden encontrarse en Chaudhry (1987), Wylie 
y Streeter (1983), y Torres, Besancon y Georges 
(2008). En este trabajo se emplea el Método 
de Diferencias Finitas por su simplicidad 
y porque la estructura de los modelos 
resultantes de dimensión finita son sistemas 


z=L del ducto. Además, se necesitan definir las 
condiciones iniciales que expresan los perfiles 
espaciales de flujo y la presión al instante inicial 
(t£ =0), esto es: 


A(2,0)=H* (2),0(2,0) =Q*(2) (6) 


Como condiciones de frontera aguas arriba 
y abajo, respectivamente (z = 0) y (z = L), se 
utilizan condiciones de Dirichlet expresadas 
como: 1) presión aguas arriba, H(0,t) = H,(t); 
2) presión aguas abajo, H(L,t) = H.(t); 3) flujo 
aguas arriba, Q(0,t) = Q.(t); 4) flujo aguas abajo: 
QU) =Q.0. 

Existen diferentes enfoques de diferencias 
finitas (ver, por ejemplo, Leveque, 2007). Así 
pues, a continuación se aproximan las deri- 
vadas espaciales parciales por diferencias 
finitas de primer orden: 


0H (2,t) _ AH; (+) _ H,(t)-H,,1(t) 


Hessenberg (Bernard, Sallet, $ Sciandra, 1998), 9 Az, Az, 
característica deseada para estimar los estados 0Q(2 t) AQ; (t) Q;(+) a (+) 
del fluido H(z,t) y Q(z,t), con observadores de de ES ma = dE - 


alta ganancia. 
Modelos finitos 


Para obtener un modelo de dimensión finita a 
partir del conjunto (4) y (5), es necesario definir 


1 1 


y se definen las derivadas temporales de la 
siguiente forma: 


condiciones de frontera que expresen los perfiles TN Ñá 
temporales de H(z,t) y Q(z,t) en los bordes del 9Q (2, t) 90, (+) 
ducto, i. e. en las coordenadas espaciales z = 0 y al :=Q; 
+ Cuadro 1. Parámetros y variables físicas del modelo del fluido. 
2 
a 
E Símbolo Parámetro/variable Unidad 
So 2 Gravedad (m/s?) 
Z b Velocidad de la onda de presión (m/s) 
+ $ Diámetro (m) 
E JE Longitud (m) 
> A, Área (m?) 
_ y Coeficiente de fricción - 
> Q (71) Flujo (m/s) 
H (z,t) Carga de presión (m) 
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donde el índice ¡ denota la sección espacial 
discretizada. Para reducir la presentación 
del modelo se omite la variable + en casos 
innecesarios y se usa la notación: 


H(z,t) = H,(t):= H, presión aguas arriba 
de la sección i 
Q(2,t) =Q;(t):=Q; flujo en la sección ¡ 


Así, en cada sección i del ducto, el conjunto 
de ecuaciones del fluido se transforma en: 


Fl, =a, LiRia (7) 
Az; 
so HH, 
rd a (5) 


1 


con 4,=84, a, E 
ZA, 294, 
Por otro lado, para considerar en las 
ecuaciones (7) y (8) el efecto de una fuga en la 
posición z ES debe incluir el flujo de extracción: 


Qy(2,,t)=H,()0; H(z,1t) (9) 


donde of = /234/C; =0;Ar es el área del 
orificio; Cr el coeficiente de descarga; Hp la 
función Heaviside, que representa la ocurrencia 
de una fuga al instante £,. 

Para el caso de varias fugas, se define n, 
como la cantidad total de fugas en una tubería; 
a 11, COMO el identificador de cada fuga, 
y Q,(2,,t) como el flujo de fuga asociado con 
k. Así, la presencia de la fuga k al instante ta 
se describe por H,(0, donde el coeficiente 
O, Y su posición Zs, caracterizan cada evento. 
Cuando se presenta una extracción Q ¡ entre las 
secciones Az, y Az, ,, la conservación de la masa 
se expresa como Q,- Q,,, — Q, = 0. Así pues, 
en cada punto de fuga, la relación de los flujos 
entre secciones: 


¡+1 


AQ; _ 


Q; E Qi 
Az. 


1 1 


Az 


se reemplaza por: 


AQ, _Q-Qn-0, 
Az, Az, 


1 1 


con Qs expresado por (9). De esta manera, el 
ducto con extracciones se describe como: 


q 
Q; A A Qi CO) 


Ha 0-04) (11) 


para ¡=1.,..., n y con variables exógenas H, = H, 
y H,,., = H, que son las condiciones de frontera. 


n+1 


Sistema de diagnóstico para el ducto 


A partir del modelo analítico del fluido (10) y 
(11), con Q, = 0, y considerando como únicas 
mediciones disponibles los flujos y las cargas 
de presión en los extremos del ducto, es decir, 
H, H.y Q, y Q,, a continuación se presenta 
el desarrollo de los algoritmos de diagnóstico 
para cinco escenarios anormales: 


e El: falla en la bomba aguas arriba. 

e  E?2: falla en la bomba aguas abajo o en el 
sistema de almacenamiento. 

e  E3: falla en el sensor de flujo aguas arriba. 

e  E4: falla en el sensor de flujo aguas abajo. 

e  E5: fuga en una posición desconocida del 
ducto. 


La figura 3 muestra el diagnosticador de 
fallas propuesto y formado por dos etapas: la de 
detección y aislamiento, y la de reconstrucción. 
La primera consta de cinco relaciones 
redundantes para aislar las fallas y la segunda 
etapa, diseñada vía observadores dinámicos, 
reconstruye la falla una vez identificada. 
Cabe señalar que sólo un observador se activa 
dependiendo del síntoma identificado. 


Detección y aislamiento de fallas 
Para obtener residuos ante los escenarios de 


fallas (E1-E5), se parte del modelo nominal (10) 
y (11), considerando las fallas como elementos 
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Detección y aislamiento de fallas 


Sirl,1r2, 14,150 El 
Si rl, 13, r4 y 15 0=E2 
Sir2, 13, 14 0 =E3 
Sir2, 13, 15 *0=)E4 
Si rl, 12, 13, 14,150 )E5 


Observadores para reconstrucción de fallas 
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Figura 3. Esquema de diagnóstico de fallas para un ducto. 


aditivos tanto en sensores como en bombas y 1 =0,-Q, (12) 
fugas. 

El residuo más usado es muy simple y se el cual es cero si el ducto opera en condiciones 
genera con la diferencia de flujos: normales y se desvía de cero cuando existen 
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fugas en el ducto o alguna de las mediciones 
del flujo es errónea; es decir, r, se desvía de cero 
ante los escenarios E3, E4 y E5. 

Considerando el ducto de longitud L en 
estado permanente, a partir de (10) y (11), la 
relación: 

1, =uQ, [Q.++(H. -H,) (13) 
vale cero en condiciones normales y se desvía 
de cero con cuatro escenarios de falla: El, E2, 
E3 y E5. 

Dado que (13) no depende de la medición 
del flujo aguas abajo Q, ésta es insensible 
a mediciones en dicha variable. De manera 
análoga, por simetría del modelo, el residuo: 

=10,[0,)+2(H.-H,) (14) 
vale cero cuando el ducto está libre de fugas 
y la información de la presiones y del flujo 
aguas abajo son correctas. Por el contrario, su 
valor se desvía de cero cuando existe una fuga 
o alguna de las mediciones mencionadas es 
errónea. Dado que r, no depende de Q, éste es 
insensible a mediciones en el flujo de entrada. 

Con ayuda de los residuos (12-14) se 
obtienen las tres primeras filas de la matriz 
de fallas (17), donde una desviación en una 
variable es antecedida de O, y * denota un 
residuo diferente de cero con la anormalidad 
correspondiente en la columna. Se observa 
de esta matriz que si se asume solamente un 
escenario presente, es posible detectar qué 
sensor de flujo está dañado tanto en condición 
de fuga como normal. Por el contrario, errores 
en H, y H, no pueden aislarse, dado que ambos 
residuos, r, y r,, son sensibles a ambos errores. 

Considerando el modelo discretizado con 
dos secciones de igual tamaño se obtiene el 
sistema de ecuaciones: 


del que eliminando la variable desconocida H , 


se tiene la ecuación integro-diferencial: 


1, =Q,-HQ, [Q.|+0, $, (Q.-Q,)dr-0,H, (15) 


_ AA, 
Con 21 y a - . Así pues, r, vale 
cero ante Ea en la presión aguas 
abajo, H, y fugas en el cuerpo del ducto. Por 
simetría del modelo se tiene la ecuación del 
residuo: 
t 

15= Q, -1Q, ¡AAN (en -Q,)dt-a,H, (16) 
que no depende de la presión de entrada H.. 
Anexando los dos nuevos residuos (15) y (16) a 
la matriz de firmas de casos anormales se tiene: 


óQ, SQ, 9H, SH, Q, 


n .. . 
ro. . . 0. 
ñ . «. 0 
.. . . . 
... . . (17) 


En esta matriz se observa que todos los 
residuos son sensibles ante la presencia de 
la fuga (última columna) y teóricamente 
cualquiera de ellos puede ser usado para su 
diagnóstico. Sin embargo, debido a que no 
todos los residuos responden de la misma 
forma ante la fuga, se recomienda usar el 
residuo cuya magnitud sea más sensible. 
Para las mediciones de los gastos y señales de 
acción asociadas con las presiones se tiene un 
patrón de residuos distinto para cada escenario 
de condición anormal, por lo que las fallas 
consideradas son todas aislables. 


Algoritmos de reconstrucción de fallas 


Como segunda etapa para reconstruir las cinco 
fallas propiamente dichas, se parte del modelo 
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no lineal de la falla correspondiente, en donde 
los parámetros 0(t) asociados con las fallas son 
estados adicionales de acuerdo con el vector 
(3) y pueden estimarse vía algoritmos de 
observación no lineales. En este trabajo se usa 
el algoritmo de observación de altas ganancias 
descrito en el Apéndice A. 

El punto de partida del algoritmo con el 
modelo discretizado del fluido (10) y (11) 
consiste en llevarlo a la forma: 


Fa(1):u(0,0(9)) =F(x.(0)2)=Gv (9) 


Obteniéndose: 


3,(t) =F(x,(t),u(+))+G(y(+)) (18) 
donde x (t), el vector de estados extendido, está 
definido por el estado x(t) y los parámetros 0(t) 
a estimar y que están asociados con la falla en 
cuestión, mientras que el vector de entradas 
está dado por: 


u(2) [un (9)-22(0)] -[H.(9)-4,(0]'. 


conformado por las presiones en los extremos 
de la tubería. En el caso de operar el ducto 
con presiones constantes o dinámica lenta, se 
asume u(t) =1,. 

Para el diseño de los observadores no 
lineales, se considera que las variables que 
pueden medirse son los flujos en los extremos 
de la tubería. Es decir: 

T 
y=h(x,(1))=[0.0,] (19) 

Sin embargo, no en todos los casos de 
reconstrucción de fallas se requieren ambas 
mediciones, así que para cada escenario de falla 
se definen las variables medidas utilizadas. 


Reconstrucción de una fuga en el ducto 


Asumiendo que la instrumentación está en 
buen estado, el objetivo ahora es estimar la 


posición y el coeficiente de una fuga en el 
ducto. A partir de (10) y (11), en presencia de 
una fuga, el modelo debe dividirse en al menos 
dos secciones con diferentes flujos: el previo a 
la fuga y el posterior. Por esta razón, se necesita 
un modelo de dimensión finita del fluido que 
represente al menos ambas secciones. 

Se propone entonces el sistema de ecua- 
ciones: 


O, =Q, IQ, l+2(8,-81,) 


H, e Q, -0,./H,) 


Ó =10 lQ (8, -H,) 


2 


(20) 


donde la posición de la fuga está dada por 
el tamaño de la primera sección Az, E [0, L], 
mientras que la segunda sección es de tamaño 
Az, = L-Az.. 

Para completar el modelo, al vector 
de estados del sistema (20) por ejemplo 
x(+)-[0, H; Q.). se le agregan dos estados 
adicionales asociados con los parámetros de 
una extracción desconocida o fuga: 


o(+)=[az, 0,] 


formando así el vector aumentado (3). Tomando 
en cuenta que estos parámetros son constantes, 
sus derivadas se igualan a cero, por tanto: 


Áz, =0, 9 =0 


Siguiendo el procedimiento del observador 
de alta ganancia descrito en el Apéndice A, 
para el vector de estados: 


x.(+)=[0, HQ, Ar; as] 


los parámetros asociados con la fuga pueden 
ser reconstruidos con la estimación de 0(t), 
utilizando el observador: 
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A A a; A 
o | "0.10 l+ (5.4) 
a 2 [6 00,8 
E AN e ls Uf 
Íñ, Áz, il 
A Ñ A UN A 
es => 10,19,1+7 7 -(8,-5,) 
e 0 
O; 0 

-K(.) Q.=Q, (21) 
Q,-Q 


22 0 
MÁZ 
En A, 0 
K(.)= 0 3h 
3,(L-áz,) 
a 
0 rn | 
mn. 3a 3 3, (L-áz,) 
H. 2a,H, 
Mz, (L-áz,) 
4,4, ñ, 


donde A es un parámetro de diseño para ajustar 
el tiempo de convergencia de los valores 
estimados a los valores reales (ver Apéndice). 
Este observador es concebido considerando 
como vector de salidas (19), es decir, las 
mediciones de los flujos en ambos extremos. 


Reconstrucción de fallas en los medidores de flujo 


Para los escenarios de fallas en los medidores 
de flujo, se supone que éstas son aditivas del 


tipo sesgo, conocido como offsets en inglés, así 
que los modelos de falla se reducen a: 


Q. =Q,+00, 


Q, =Q,+80, 

Sustituyendo estos modelos en (10) y (11) 
para una sola sección de longitud L se obtienen 
dos modelos extendidos: 

e Para fallas en el sensor aguas arriba: 

O. =1.(0, -9Q,)IÓ, -5Q, |l+2(A, -H,) 

30, =0 
conzo(1)=[6,50,]' 

e Para fallas en el sensor aguas abajo: 


O, =1u(0, -90,)/0, -90, 


Lot qa 
+ pH, H,) 
30, =0 


T 


con x.(t)=[0,5Q,] 


En esta pareja de modelos, las variables des- 
conocidas 30, =0 y JO) =0, respectivamente, 
representan la dinámica constante del sesgo 
a estimar. La elección de un modelo finito 
discretizado con una sola sección obedece a 
que a lo largo del ducto el flujo es el mismo en 
estado permanente en ausencia de extracciones. 

Siguiendo el procedimiento descrito en 
el Apéndice A, la ecuación del algoritmo de 
observación de alta ganancia diseñado para el 
caso de falla en el medidor de flujo aguas arriba 
se reduce a: 


-KC|O -0,] (22) 
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con la ganancia: 


24 
K()=| -22 
() LAN 
Q 
2110, -50, +10,sign[0, -50, 


190, sign [9 - 50) 


Suponiendo que (6, -90,) > 0, la ganancia 
se reduce a: 


2h 
Ka E 


2/0, -50,) 2 


En este caso, el diseño del observador utili- 
za solamente la medición del flujo aguas arriba, 


y = h(x.(t)) = Q.. 
Reconstrucción de fallas en sistema de bombeo 


En el caso de fallas en los sistemas de bombeo o 
almacenamiento en los extremos del ducto, los 
modelos de falla que se consideran son: 


Á, =H,+0H, 
á,=H,+0H, 


Para diseñar los observadores que estimen 
las fallas de presión en los extremos se necesita 
un modelo discretizado al menos con dos 
secciones espaciales de manera similar que para 
la reconstrucción de una fuga. Para el caso de 
fallas aguas arriba, el modelo que se considera 
para el diseño del algoritmo de observación es: 


O, =1Q, 


Q 


a po 
rea -dH, -H,) 
H,--E-(Q,-0,) 
f Az, e Ss 


O, =1Q, 


Q.|+ 


a 
2 (4,4) 


SH, =0 


Por simetría, para una falla en la presión 
aguas abajo, el observador se basa en el modelo: 


O, = 0, lO, 


+¿¿(H,-H,) 
1, -H0.-Q,) 


Q, =1Q. lO, 


a - 
+2/H,-H,+0H, 
Az, f ] ) 
sH,=0 
El algoritmo de observación de alta ganancia 


diseñado para el escenario de falla en el sistema 
de bombeo aguas arriba está dado por: 


ó, 7 10.0,+2-(A.-A,-0H,) 
A, |. =(9,-0) 
o, 30 (a,-4 
sá, ASA! e 
j 0 
-K() 272. (23) 
Q,-Q 


el cual es concebido considerando como vector 
de salidas a las mediciones de los flujos en 
ambos extremos, por ejemplo, utilizando la 
ecuación (19). Suponiendo que el flujo satisface 
la condición Ó. > 0, la ganancia reducida del 
observador se expresa como: 


2% 0 ] 
o =N 4u0 Az,A 
a a 
K(.)= ds 
En e Az,N ¿HMQ,Az,A 
a a 


AZ ñ 410 Az,» 
a a 


O= 
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Modelo físico piloto 


La instalación experimental para detección 
de fugas es una tubería de acero de cuatro 
pulgadas de diámetro, cédula 40, en forma de 
espiral sobre un plano vertical, como lo muestra 
la figura 4. La entrada del flujo está en la parte 
inferior y su salida en la parte superior. Se tiene 
una cisterna de 10 000 litros de agua potable, 
una bomba hidráulica de 7.5 HP y válvulas de 
paso en la entrada y en la salida. Los sensores 
de flujo y presión se ubican en ambos extremos 
del ducto. Aproximadamente en la parte media 
de cada tubo horizontal se haya una válvula de 
l in de diámetro para simular fugas. 

La tubería tiene una longitud física cuando 
se hacen mediciones directas sobre ella punto a 
punto, pero la presencia de codos a lo largo del 
ducto genera una caída de presión mayor a la 
correspondiente a su longitud real, por ello es 
necesario determinar la longitud equivalente 
si se supone la instalación como una tubería 


Figura 4. Vista lateral de la instalación piloto. 


recta; de esta manera se tiene una equivalencia 
entre la tubería en espiral de 169.43 m y una 
recta de 200.17 m. 

La tubería también cuenta con seis válvulas 
(VL, V2, V3, V4, V5, V6) utilizadas para emular 
fugas, dos de ellas servo-válvulas, para realizar 
experimentos repetitivos, y localizadas a 
11.535, 49.825, 80.355, 118.365, 148.925 y 186.945 
m, respectivamente, a partir del punto de 
alimentación. Los valores para los parámetros 
de la instalación se reportan en el cuadro 2. 

La instalación cuenta con tres opciones 
para medir el flujo con tres tipos de sensor: de 
propela, ultrasónico y fuerza de Coriolis. Estos 
sensores se ubican al inicio y final del ducto. 


Sistema de supervisión, control y adquisición 
de datos SCADA 


Este sistema está compuesto básicamente por 
dos elementos: uno de hardware y otro de soft- 
ware. El primero está compuesto por una com- 
putadora personal conectada a un módulo de 
adquisición de datos vía protocolo "TCP, éste, a 
su vez, está conectado físicamente a los senso- 
res, servo-válvulas y actuador de la bomba; el 
segundo es el programa de aplicación desarro- 
llado en LabVIEW?. 

El hardware: la computadora personal 
trabaja en ambiente de Windows y en ella se 
tienen utilerías para poder conectar el software 
de LabVIEWO con un sistema modular de 
adquisición de datos con protocolo Modbus 
(BK9000 de Beckhoff ) vía una red local Ethernet. 
La comunicación entre la aplicación y el sistema 
modular se hace por un servidor OPC, que es el 
intermediario entre el programa y los registros 
de entrada-salida vía comunicación TCP. Los 


Cuadro 2. Parámetros del ducto de pruebas (modelo físico). 


Símbolo Valor Unidad 
b 1 497 (m/s) 
$ 0.1047 (m) 
1 200.16 (m) 
f 2.785x 107? - 
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módulos de adquisición usan un sistema de 
4-terminales para alambre, con ello todos los 
sensores y actuadores con diferentes tipos de 
señales puedan ser conectados directamente al 
módulo correspondiente. 

El software: es un sistema de adquisición y 
monitoreo del ducto integrado con base en la 
plataforma de la compañía National Instru- 
ments? llamada LabVIEW, que es un sistema de 
programación de propósito general compuesto 
por librerías de funciones y herramientas de 
desarrollo diseñadas específicamente para ad- 
quisición de datos y control de instrumentos. 
Los programas de LabVIEW” son llamados 
instrumentos virtuales, porque en su aparien- 
cia y operación pueden imitar instrumentos en 
boga. Funcionalmente, el sistema desarrollado 
realiza actividades de lectura de la información 
de los sensores de flujo y presión, el manejo de 
la abertura de las servo-válvulas y regulación 


de la velocidad del motor de la bomba. Depen- 
diendo de la aplicación, se tiene una opción de 
sólo registrar los datos de un experimento para 
su posterior procesamiento fuera de línea o la 
opción de la detección de una fuga en tiempo 
real con alguno de los métodos de detección 
desarrollados en el Instituto de Ingeniería, 
como es el caso de los algoritmos validados en 
esta contribución. 


Importancia de los sensores en la estimación 
de la posición de la fuga 


La ubicación precisa de una fuga no sólo 
depende del método de detección y ubicación, 
sino también del tipo de instrumentos de 
medición. La figura 5 muestra tres estimaciones 
de la posición de la fuga ubicada a 49.8 m para 
una condición abierta de válvula del 100%. 
Se usó el mismo algoritmo de estimación, 


Localización de fuga (posición real = 49.8 m) 
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Figura 5. Desviación en el estimador de la posición de la fuga, según el tipo de sensor de flujo. 
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pero con distintos sensores de flujo. Se puede 
observar que además de registrarse diferentes 
ubicaciones, las desviaciones estándar son 
notoriamente diferentes, provocando incerti- 
dumbres en la ubicación de la fuga. El sensor de 
fuerza de Coriolis genera un mejor desempeño. 


Resultados experimentales 


Con el propósito de evaluar el desempeño del 
sistema de diagnóstico de fallas en ductos, se 
realizaron varios experimentos provocando 
tanto una extracción como alteraciones en las 
lecturas registradas de los sensores de flujo y 
presión. 

El evento de la extracción desconocida o 
fuga se originó abriendo la válvula V4 a los 150 
s. El efecto de la fuga puede observarse en el 
comportamiento de los flujos mostrados en la 
figura 6, ya que después del evento se desvían 
del punto de operación, a pesar de que la presión 
parece permanecer constante, ya que el flujo de 
fuga es menor al 8% y sus efectos en la carga 
son imperceptibles. En esta figura se muestran 
también las señales medidas de las presiones 
en los extremos de la tubería utilizadas por los 


observadores para realizar las estimaciones. En 
las gráficas se han superpuesto a las mediciones 
reales valores filtrados de éstas, con la finalidad 
de tener mejores resultados en la estimación. El 
filtrado se llevó a cabo utilizando un filtro pasa 
bajas de primer orden, con una frecuencia de 
corte w,= 0.01 rad /s y ganancia unitaria. 


Aislamiento de las fallas 


Considerando los datos reales del ducto piloto, 
las figuras 7 y 8 muestran la evolución de tres 
residuos Cr YT, 1) ante escenarios anormales. La 
figura 7 corresponde al caso de una extracción 
iniciada a los 100 s, en donde se comprueba 
que los tres residuos son distintos de cero una 
vez que la extracción está presente, como lo 
pronostica la matriz de firma de fallas (17). La 
figura 8 muestra dos escenarios distintos. La 
parte superior corresponde a la respuesta de los 
residuos ante una desviación en el medidor de 
flujo aguas abajo y la inferior ante un error en 
el sistema de bombeo aguas arriba. En ambos 
casos, los dos conjuntos de residuos responden 
también de acuerdo con la matriz de falla (17) y, 
por lo tanto, con condiciones lógicas, es posible 
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Figura 6. Señales medidas de flujos y presiones en los extremos del ducto con extracción a partir de 100 s. 
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Figura 7. Evolución de los residuos 1,(t), r,(t) y r,(t), con una extracción provocada a los 100 s. 


Escenario de falla en el medidor de flujo aguas abajo con un sesgo de 50% 
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Escenario de falla en el sistema de bombeo aguas arriba con sesgo de 40% 


Figura 8. Evolución de los residuos r,(t), r,(t) y r,(t), ante dos escenarios de fallas. 
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aislar la falla presente. Esto siempre y cuando e El: reconstrucción de la falla en la bomba 
exista solamente una falla. con el algoritmo (23). 
Reconstrucción de las fallas Para la evaluación del desempeño del 
localizador con datos experimentales se uti- 
Según la naturaleza de una falla, una vez que lizan las mediciones de flujo y presiones 
ésta ha sido detectada y aislada usando las ilustradas en la figura 6. 
condiciones lógicas dadas en la matriz (17), 
la reconstrucción y el valor en el tiempo de Estimación de una fuga 
la falla se determina utilizando el algoritmo 
de observación correspondiente de acuerdo El observador para la estimación de la fuga 
con el escenario identificado. A continuación (21) se puso en marcha al tiempo t = 160 s, una 
se presentan los resultados experimentales vez que el generador de residuos identificó 
de los observadores de reconstrucción con los que el síntoma se debía a una extracción 
siguientes tres escenarios de fallas: desconocida. Cabe decir que el observador 
sólo opera de manera adecuada en presencia 

+ E5: reconstrucción de la fuga f, con el de fuga, de lo contrario no se puede garantizar 

algoritmo (21). la convergencia del error (Torres, Verde, 
e E3: reconstrucción del sesgo %Q, en el Besancon, € González, 2012). 

medidor de flujo aguas arriba con el La figura 9 muestra la evolución del estima- 

algoritmo (22). dor de la posición realizada por el observador 

mio! MIME MONIAINN JUNE TUNA PAGINE NIDTIDE: TEE TIE: TINE TS E CE E TI E E AE 
: AÑ Ñ 3 ¿ ¿ J : j ¿ ¿ ¿ ¿ 7 
so AY JO A PA 


: é 3 A j : 1 :B | | Y 
A 0% +» a E 0 IAS € a 0 A AS 


(m) 


: q : : : Posición real de fuga! E j : : : 
A | STE ES EIN OS EN AS bario AN AS bno pa 


HEEE 


50 mi | 1 | mi | L | 1 j mi | 1 | L 1 | 1 | ñ j mi j L | L 


Figura 9. Estimación con diversos valores de A. 
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Figura 10. Estimación de la posición con/sin filtro. 


para diferentes valores de A, observando que el 
tiempo de convergencia disminuye si el valor 
de A se incrementa. Sin embargo, el ruido de 
medición también aumenta con el incremento 
de A. En la figura 10 se observa la estimación 
de la posición de la fuga cuando el parámetro 
de convergencia se fija A =3 para dos diferentes 
casos, con señales filtradas y sin filtrar. En 
ambas figuras se observa que los valores 
estimados están cercanos al valor real de la 


posición real de la fuga. 
Estimación de sesgo en un medidor de flujo 


Para evaluar el observador de reconstrucción 
de fallas en el sensor de flujo aguas arriba (22), 
se sumó una señal constante a la medición 
aguas arriba de 3 x 10?m?/s. El observador fue 
inicializado en el tiempo t=0 s del experimento 
mostrado en la figura 6, es decir, cuando 
no hay presencia de fuga. Las condiciones 
iniciales del observador se eligieron diferentes 


a las de la tubería real y se sintonizó con A = 2. 
Las señales de entrada se filtraron con dos 
diferentes frecuencias de corte, cuyos valores 
se muestran en la figura 11, donde también se 
puede apreciar el resultado de la estimación. 
Es notable que si no se utilizan filtros en las 
entradas, la estimación del parámetro hereda 
el ruido de éstas. También es interesante que 
la frecuencia de corte elegida para los filtros 
incide en la estimación, porlo que hay que tener 
en cuenta este factor durante la calibración del 
estimador. 


Estimación de fallas en la bomba 


Para la evaluación del observador que estima 
fallas en el sistema de bombeo (23) se sumó 
una señal constante a la medición de la pre- 
sión aguas arriba de 1.6x 10 m. El observador 
se inicializa cuando no hay extracciones. Las 
condiciones iniciales del observador se eligie- 
ron diferentes a las de la tubería y se sintonizó 
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Figura 11. Estimación de offset en el medidor de flujo aguas arriba. 


con A = 1. En la figura 12 se puede apreciar el 
resultado de la buena estimación realizada por 
el observador con y sin filtros. 


Conclusiones 


En este artículo se presentó un sistema robus- 
to de diagnóstico para fallas en ductos imple- 
mentado en un modelo piloto del Instituto de 
Ingeniería de la UNAM. Este sistema actual- 
mente está programado en LabVIEW” para la 
adquisición de datos y en MATLAB? para el 
diagnóstico y la reconstrucción de fallas. La 
originalidad de este sistema es su capacidad de 
detectar, aislar y reconstruir diferentes tipos de 
falla en ductos, novedad y aportación principal 
del presente trabajo. 

La principal característica de este sistema 
de supervisión es que está compuesto de dos 
etapas con tareas específicas. En la primera 
etapa, el objetivo es distinguir entre cinco 
escenarios de fallas, entre los que se encuentra la 
fuga del ducto. Esta tarea es realizada mediante 
el uso de relaciones redundantes obtenidas a 


partir del modelo de transporte del fluido en 
el ducto. En la segunda etapa, la tarea es la 
reconstrucción de la evolución en el tiempo de 
la falla que se ha identificado. Los resultados 
de la evaluación del sistema de diagnóstico 
con datos experimentales fueron satisfactorios, 
haciéndolo viable para su implementación en 
ductos reales. 


Apéndice 
Diseño de un observador de alta ganancia 


Gauthier et al. (1992) demostraron que a cual- 
quier sistema representado por (1), que posea 
la propiedad de poder ser transformado en la 
forma triangular siguiente: 


(24) 


se le puede diseñar un observador de alta 
ganancia con: 
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Figura 12. Estimación del error en el sensor de presión aguas arriba. 


o 1 0 
0.0 1 
A= 0 
: 1 
0 : 0 


(25) 


El procedimiento de diseño para un 
observador, considerando el caso particular del 
modelo en espacio de estados (18) con entradas 
constantes 1, que expresa de manera compacta 
el modelo finito de transporte de fluido en un 
ducto dado por las ecuaciones (10) y (11) se 
reduce a los siguientes pasos. 


Paso 1: aplicar la transformación no lineal 
D(x,(t)) = E(t) al sistema (18) mediante el 
siguiente cambio de coordenadas: 


donde L, es la derivada de Lie de la salida h(x,) 
t)) con respecto al operador f (consultar Isidori, 
1995, para profundizar). 

Paso 2: construir la estructura del algoritmo 
de observación: 


E(t) - AE(1)+p(8(8),u(1) 


Sistema transformado 


s”C*(CE(+)- y(+)) 


Término de corrección 


(26) 


que no es más que una copia del modelo 
del sistema transformado, y un término de 
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corrección que asegura la convergencia del 
algoritmo con la matriz S, el parámetro de 
diseño. 

Paso 3: regresar el observador (26) a las 
coordenadas originales, obteniendo: 


incluye la matriz S, que es simétrica, definida 
positiva, que asegura la convergencia del 
observador, siendo la solución única de la 
siguiente ecuación matricial: 

=AS+ A'S+SA-C'C=0 (28) 
donde A > 0 es un parámetro de diseño; las 
matrices A y C están dadas en el modelo 
canónico (24). 

Nota: Gauthier et al. (1992) demostraron 
que el error de estimación converge exponen- 
cialmente a cero si el parámetro A se elige sufi- 
cientemente grande. Además, la sintonización 
del observador se reduce a la calibración de 
este simple parámetro y la elección de un valor 
de A grande asegura una convergencia rápida 
de los estimados a los valores reales, de aquí 
el nombre alta ganancia. Sin embargo, una ga- 
nancia muy alta produce sobrepasos durante 
la convergencia y amplificación del ruido. Por 
ende, la elección de la ganancia debe conside- 
rar estos factores. 
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